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^ ■ The photoproduction of large-pT charged hadrons and of prompt pho- 

O ' tons is discussed, for the inclusive case and with an associated jet, using 

, predictions from the NLO partonic Monte Carlo program EPHOX. Com- 

. ■ parisons to recent HERA data are also shown. 

^ ■ 1 Introduction 

o ■ 

. At HERA, the photoproduction of jets has been measured with high statistics 

' over the past years ^ i2, ^ . Data on and single charged hadron 4 and prompt 

. photon [SI E] production are also available and have been compared to theo- 

retical predictions 115| . where in the prompt photon case the statistics is of 
course lower as the cross sections are small. 

In photoproduction reactions, an almost real photon, emitted at small angle 
from the electron, interacts with a parton from the proton. The photon can 
either participate directly in the hard scattering ("direct photon") or act as 
Oh. a " resolved photon" , in which case a parton stemming from the photon takes 

part in the hard interaction. Therefore photoproduction reactions offer a unique 
opportunity to constrain the parton distributions in the photon, especially the 
gluon distribution g'^{x). The latter is rather poorly constrained by LEP 7*7 
data as it enters in this reaction only at next-to-leading order (NLO), whereas 
^ ' in ep photoproduction it enters already at leading order. 

On the other hand, photoproduction reactions at HERA also could serve 
to constrain the gluon distribution in the proton. To determine the latter to a 
better accuracy is of particular interest in view of the LHC with its large gluon 
luminosity, gg — > H being the dominant production mode for a light Higgs 
boson. As the photon represents a rather "clean" initial state, photoproduction 
reactions can be used to probe the proton in a way which is complementary to 
other experiments. 

In what follows, the kinematic regions where gluon initiated subprocesses 
play a significant role will be identified and it will be argued that their relative 
contribution to the cross section can be enhanced by appropriate cuts, for the 
proton as well as for the photon case. For a more detailed study, the reader is 
referred to JSI ■ 

Comparisons to HERA data, especially to very recent (preliminary) HI data 
on prompt photon plus jet production, will also be shown. 

^ To appear in the proceedings of the Ringberg Workshop " New Trends in HERA Physics 
2003", Ringberg Castle, Tegernsee, Germany, September 28 - October 3, 2003. 
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2 Theoretical aspects of charged hadron and prompt 
photon production 

We will concentrate here on the two reactions "fp (+jet) +X and ■yp 

7 (+ jet) + X. Identifying a jet in addition to the prompt photon or hadron allows 
for a more detailed study of the underlying parton dynamics, and in particular 
for the definition of observables Xobs, xll which are suitable to study the parton 
distribution functions. 

Comparing the reaction of charged hadron production versus prompt pho- 
ton production, one is tempted to say that prompt photon production is more 
advantageous: Due to photon isolation, the dependence of the theoretical pre- 
diction on the non-perturbative fragmentation functions is negligible. Further, 
as will be shown in the following, the NLO prediction for the prompt photon 
cross section is not very sensitive to scale changes, which cannot be said for 
the {+ jet) cross section 0. On the other hand, the 7-)- jet cross section is 
orders of magnitude smaller than the + jet cross section, and the identifica- 
tion of the prompt photon events among the huge background from the decay 
of neutral pions is experimentally not an easy task. 

Photon isolation 

Prompt photons can originate from two mechanisms: Either they stem directly 
from the hard interaction, or they are produced by the fragmentation of a hard 
parton. In order to distinguish them from the background stemming from the 
decay of light mesons, an isolation criterion has to be imposed: The amount 
of hadronic transverse energy E^^'^ , deposited inside a cone with aperture R 
centered around the photon direction in the rapidity and aziniuthal angle plane, 
must be smaller than some value E^^^: 



For the numerical results we will follow the HERA conventions E^^^ = 0.1 p^ 
and R = 1. 

Apart from reducing the background from secondary photons, isolation also 
substantially reduces the fragmentation component, such that the total cross 
section depends very little on the photon fragmentation functions. 

Factorisation and scale dependence 

The cross section for prompt photon production can be written symbolically as 
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In the case of the production of a single hadron, the partonic cross section 
da'^" is of course zero, and the fragmentation functions D^/c{z, Mp) have to be 
replaced by the ones for hadrons, Di^/^{z, Mp). 

In photoproduction, _F'^/g(a;e, M) is a convolution 



F,/e(xe,M) = j dx^ j dy5{x'^y-Xe)f^/e{y)Fa/^{x\M) 



(2) 



where the spectrum of quasi-real photons emitted from the electron is described 
by the Weizsacker- Williams approximation 
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The function Fa/~f{x'' , M) in eq. 10) de- 
notes the parton distribution function 
for a parton of type "a" in the re- 
solved photon. In the case of a direct 
initial photon, one has Fa/yix"^ , M) = 



scale dependence of a'°' for 7 + X 



Therefore one can try 



to switch on/off the resolved photon by 
suppressing/enhancing large x'^ . 

As the perturbative series in as 
for the partonic cross section is trun- 
cated, the theoretical prediction is scale 
dependent, depending on the renor- 
malisation scale fi as well as on the 
initial/final state factorisation scales 
M /Mp. While the leading order cross 
section depends very strongly on the 
scales, the NLO cross section is already 
much more stable, as can be seen for 
example from Fig. ^ for the prompt photon inclusive cross section. 



Figure 1: The scales have been set to 
fjL — M — Mp = C p7p with C varied 
between 0.5 and 2. 



Observables x^^x^ 

In order to reconstruct the longitudinal momentum fraction of the parton stem- 
ming from the proton respectively the photon from measured quantities, one 
can define the observables 



^obs 

r-P 



PT e-" + -Bf ' 



2Ef 
Pt + -By * e' 
2EP 



(3) 



where pt is the momentum of the prompt 7 or the hadron, and 77 its (pseudo- 
)rapidity. However, as the measurement of E^r^*" can be a source of systematic 
errors at low Ep values, the following variable, which does not depend on iJ^*^*, 
might be more convenient: 



Pt (e 



2i?P.T 



(4) 



3 



The variable also has the advantage that it has a smoother behaviour for 
a:T 1. 

3 Numerical studies 

Unless stated otherwise, the following input for the numerical results is used: 
The center of mass energy is ^/s = 318 GcV with = 27.5 GeV and Ep = 
920 GeV. The maximal photon virtuality is Q^ax = 1 GeV^, and 0.2 < y < 0.7. 
For the parton distributions in the proton we take the MRSTOl 17" parametri- 
sation, for the photon we use AFG04 18 distribution functions and BFG [T??j 
fragmentation functions. For the charged hadron fragmentation functions we 
use BFGW(23j as default. We take nj = A flavours, and for as(/x) we use an 
exact solution of the two-loop renormalisation group equation, and not an ex- 
pansion in log(/i/A). The default scale choice is M — Alp = i^l = px- Jets are 
defined using the /c-r-algorithm. 

3.1 The NLO program EPHOX 

All numerical results shown here are obtained by the NLO partonic Monte Carlo 
program EPHOX 20 , which allows to obtain integrated cross sections as well 
as fully differential distributions for the reactions considered here. The program 
contains the full NLO corrections to all four categories of subprocesses direct- 
direct, direct-fragmentation, resolved-direct and resolved-fragmentation. It also 
contains the quark loop box diagram 717 — > 75 with a flag to turn it on or off. 
EPHOX can be obtained at http : //wwwlapp . in2p3 . f rT lapth/PHOX JAMILY /main . htinl| 
together with a comfortable user interface and detailed documentation. 

3.2 The gluon distribution in the photon 

In order to constrain the gluon distribution g'^{x'^) in the photon, one has to 
focus on a kinematic region where the gluon content of the resolved photon 
is large, i.e. on small values of x'^ . According to eqs. © and (0J, small 
corresponds to large values of rj and rj^'^^ , which means that the forward rapidity 
region is particularly interesting. Fig.|3illustrates for the prompt photon plus jet 
cross section that the gluon distribution g'^{x'^) only becomes important in the 
very forward region if > 1.5 if the rapidity of the jet is integrated over the range 
— 1 < Ty^'^* < 2.3. However, the relative importance of the gluon distribution 
{x) can be enhanced by restricting both rapidities, the one of the photon 
and the one of the jet, to the forward region. Fig.lJt- shows that with the cuts 
^7 ,^jet y ^Yi^ resolved photon component makes up for a major part of the 
cross section, but mainly consists of quarks, while for rf > 0.5, rj^'^* > 1.5, the 
gluon in the photon contributes about 40% to the total cross section, as shown 
in Fig.Eb. 
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Figure 2: Rapidity range where the gluon in the photon becomes important 
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Figure 3: The rapidity cuts j]'^ > 0.5, ry-''^* > 1.5 substantiaUy enhance the 
relative contribution of the gluon in the resolved photon 
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3.3 The gluon distribution in the proton 



hadron + jet 

-2 < v"''" < 0.5 
1 E,'" > 5 GeV, p," > 7 GeV 



total 

only 

MRST01 pdfs 
BFGW ffs 



The gluon distribution in the proton 

is large for small x^, corresponding £ 

to small rapidities. From Fig.0] one '■.^ 

can indeed see that the relative con- \ 

tribution of subprocesses initiated by % 

a gluon from the proton to the + 

jet cross section is about 73% if both 

rapidities are restricted to the back- 
ward region, —2 < r],?]^^* < 0.5. The 

analogous is not true for the prompt 

photon -I- jet cross section, where the 

gluon contribution makes up only 13% 

of the total cross section if the rapidi- 
ties are restricted to the range — 2 < 

7?'',77^'^* < 0. This is due to the fact 

that at small values, x'^ is large, 

such that direct initial photons should 

dominate. However, the subprocess 

+ 'f ^ 7(dircct) +jet does not exist at leading order, and the subprocess 
+ J ^ g-|-jet, where the quark subsequently fragments into a photon, is 

suppressed by isolation. Therefore, in the case of the 7 -|- jet cross section, the 



Figure 4: The gluon (from the proton) 
contribution to the h^+jet cross section 
at small rapidities. 




Figure 5: The relative contribution of gluon initiated processes to the 7 + jet 
cross section can be enhanced by the rapidity cuts < rp ,7]^'^* < 2. The px 
cuts are E^""* > 5 GeV, > 6 GeV. 



subprocess gP + q'^ j + q is the dominant one involving gP whereas in the case 
of the hadron+jet cross section, gP + j ^ q + q is dominant at small x^ . The 
virtue of this behaviour of the prompt photon cross section is that the region 
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where initiated subprocesses are important is not restricted to negative ra- 
pidities. As can be seen from Fig.O the rapidity cuts < rj'^'^^* < 2 enhance the 
relative contribution of the gluon gP to the 7 + jet cross section, because they 
select a region where g^ + initiated subprocesses are important. Note that 
this rapidity domain is more accessible experimentally than the very backward 
region. 



4 Comparison to HERA data 



4.1 Charged hadrons 

For the case of single charged 
hadron production, only data for 
/i^+X, but not for /i^+jet -|-X are 
available so far 0] . The analysis 
uses a minimum pT of 3 GeV for 
the hadron, which is rather close to 
the non-perturbative regime. This 
induces a large theoretical uncer- 
tainty due to differences in the 
parametrisation of the fragmenta- 
tion functions j21ll22ll^ at low pT 
(see Fig.EJ, and also a large scale 
dependence. If a minimum pT of 
7 GeV is chosen for the hadron, 
the cross section is reduced, but 
the theoretical uncertainty is much 
smaller [5] . 




Figure 6: Predictions obtained with differ- 
ent fragmentation functions compared to 
HI data. 



4.2 Prompt photons (inclusive) 

For inclusive prompt photon pro- 
duction, ZEUS data as well as 
preliminary HI data [J] are avail- 
able^. In Fig.CI one observes that 
the ZEUS data are above the NLO 
prediction in the backward region. 
At large rapidities, there is a trend 
that theory is above both HI and 
ZEUS data. As already men- 
tioned in the introduction, this ef- 
fect is very likely due to photon 
isolation: A sizeable amount of 
hadronic transverse energy in the 
isolation cone may stem from the 
underlying event. Therefore, even 



1996-2000 data 




5 GeV < p/ < 10 GeV 
0.2 < y < 0.9 



Figure 7: Comparison of inclusive 
photon data with the NLO QCD 
tion. 



prompt 
predic- 



^In order to be able compare to the ZEUS data taken at Ep = 820 GeV and 0.2 < j/ < 0.9, 
the HI data taken at Ep = 920 GeV and 0.2 < y < 0.7 have been corrected using PYTHIA Q. 
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direct photon events may be rejected, leading to a decrease of the cross 
section. This fact cannot be simulated by a partonic Monte Carlo, such that 
the NLO predictions tend to be above the data in kinematic regions where the 
underlying event activity is expected to be large. 



4.3 Prompt photon + jet 



In the rapidity distribution 
for the 7 + jet cross sec- 
tion, the effect mentioned 
above is again visible at 
large rapidities. To esti- 
mate its impact, HI re- 
cently made an analysis [7| 
where the NLO result is 
corrected for multiple inter- 
action effects. As can be 
seen from Fig.|Sl the correc- 
tions are indeed most pro- 
nounced in the forward re- 
gion where the underlying 
event activity from the re- 
solved photon remnants is 
expected to be larger, and 
they improve the agreement 
between data and theory. 
Fig.|Sl also shows that the 
7-l-jet cross section is very 
stable with respect to scale 
variations. 
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NLO + Ml 
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-1 -o.s -0.6 -0.4 -0.2 0.2 0.4 0.6 o.a 



Figure 8: Comparison of preliminary HI data to 
the NLO prediction with different scale choices. 
NLOH-MI denotes the result where the partonic 
prediction has been corrected for multiple inter- 
action effects. 



4.4 Study of intrinsic (A;^) 

The ZEUS collaboration made an analysis on prompt photon -I- jet data to 
study the "intrinsic" parton transverse momentum (fcy) in the proton 6 . To 
this aim, the (fc^) - sensitive observables p±, the photon momentum compo- 
nent perpendicular to the jet direction, and A0, the azimuthal acoUinearity 
between photon and jet, have been studied. To suppress contributions to (/ct) 
from the resolved photon, the cut a;"*"* > 0.9 has been imposed. To min- 
imise calibration uncertainties, only normalized cross sections have been consid- 
ered, as an additional (kx) mainly changes the shape rather than the absolute 
value of the cross section. The best fit with PYTHIA 6.129 lead to the result 
(fcr) = 1.69 ± 0.18(stat) ^Q2Q (sys) GeV, which includes the shower contribution 
to {kr). On the other hand, the NLO program EPHOX is able to describe 
the data without including an additional (fcr)- This is shown in Fig.^ where 
also the minimum has been varied in order to estimate the impact of an 
uncertainty in the determination of the jet energy and of choosing symmetric 
cuts {p'^ > 5 GeV has been used in the experimental analysis). 
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Figure 9: ZEUS data for the normalised cross sections differential in p± and 
A(j) compared to the EPHOX prediction for different cuts on the minimum jet 
transverse energies. 

5 Conclusions 

The reactions 7p — > 7 + jet + X and 7p — > /i^ + jet + X offer the possibility to 

constrain the gluon density in the photon and in the proton. It has been shown 
that appropriate rapidity cuts can enhance the relative contribution of gluon 
initiated processes. 

Comparing the predictions of the partonic Monte Carlo NLO program EPHOX 
to HERA data, the following observations can be made: For inclusive single 
charged hadron production, the HI data are well described, but the theoretical 
uncertainties from fragmentation functions and scale variations are large. These 
large uncertainties can mainly be attributed to the fact that a minimal pt of 
3 GeV for the hadron is too close to the non-perturbative regime. A p™™ of 
7 GeV improves the stability of the theoretical predictions. 

The rapidity distributions for inclusive prompt photon production as well as 
for 7+jet show an interesting feature: At large rapidities, the NLO prediction 
always tends to overshoot the data. It has been argued that this behaviour might 
be attributed to photon isolation: Due to hadronic activity stemming from the 
underlying event in the isolation cone, even direct photons in the final state 
may be rejected by the isolation cut, thus leading to a lower cross section than 
the (partonic) Monte Carlo prediction, especially in the forward region where 
the probability of resolved photon remnants is higher. A recent HI analysis 
which corrects the NLO prediction bin per bin for multiple interaction effects 
corroborates this explanation. 

Finally, the predictions of EPHOX have been compared to a ZEUS analysis 
on prompt photon+jet data where the intrinsic {kr) of the proton is studied, 
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and it has been found that the NLO prediction describes the data very well 
without introducing an extra (/ct)- 
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